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O . Abstract Local parametric statistics of zeros of Husimi representations of quantum eigen- 
states are introduced. It is conjectured that for a classically fully chaotic systems one 
should use the model of parametric statistics of complex roots of Gaussian random poly- 
^ ■ nomials which is exactly solvable as demonstrated below. For example, the velocities 
. (derivatives of zeros of Husimi function with respect to an external parameter) are pre- 
^! dieted to obey a universal (non-Maxwellian) distribution 

^ ' where o"^ is the mean square velocity. The conjecture is demonstrated numerically in 
a generic chaotic system with two degrees of freedom. Dynamical formulation of the 
^ ' "zero-flow" in terms of an integrable many-body dynamical system is given as well. 
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1 Introduction 

The intense research in the so-called quantum chaology has produced many different sig- 
natures of classical chaos in the corresponding quantum Hamiltonian systems. It has 
been found that energy spectra and eigenstates of classically fully chaotic quantum sys- 
tems have universal statistical properties which can be described by stochastic models 
with no free parameters, like random matrix theory (RMT). 

In the Bargmann (or Husimi) representation, eigenstates of quantum systems can be 
uniquely represented in terms of complex analytic functions in phase space variables. 
Quite recently, it has been proposed that in case of 1-dim systems (or Poincare surface 
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of section reductions of 2-dim systems |T0|)! where phase space is two dimensional, and 
one has only one complex variable z = q + ip, any eigenstate can be uniquely represented 
by the collection of (complex) zeros of its Bargmann or Husimi representation. This has 



been called the stellar representation. It has been found |10] that the structure of zeros 
of Husimi representation of an eigenstate is reminiscent of the structure of classical phase 
space: (i) in classically fully chaotic systems the zeros tend to spread uniformly over the 
whole classically allowed region of phase space [0 , (ii) in classically integrable systems the 
zeros lie on 1-dim. anti-Stokes curves (iii) while in generic mixed systems a part of 
zeros uniformly cover chaotic components of classical phase space, zeros in regular islands 
lie along 1-dim. classically invariant curves - tori, while substantial part of zeros lie on 
1-dim. classically noninvariant (generalized anti-Stokes) curves. 

In the following we shall be interested only in the case of classically fully chaotic 
systems. Let us assume that the Hamiltonian of our chaotic system can be statistically 
described by Gaussian orthogonal/ unitary (in the presence/absence of anti-unitary sym- 
metry) random matrix in a generic basis. If we choose the harmonic basis then the 
Bargmann representation of an arbitrary eigenstate is given by the entire function 



n=o V n! 



(1) 



and the Husimi representation in appropriately scaled units is 1/(2;) pe~'^'^, where z = 
q + ip. Cn are the coefficients of the expansion of that eigenstate in the harmonic basis, 
which are, by assumption and hence according to RMT, real/complex Gaussian (pseudo) 
random variables. Ideally, as RMT predicts, there should be no correlations among them 
(cnC^) = Snm, although we believe that for nonzero shortrange correlations among the 
coefficients c„ our general conclusions are still valid. Since Taylor expansion of f{z) around 
arbitrary point zq is convergent one may define and study Gaussian random polynomials 
of order N and the statistical properties of their roots and in the end, if neccesary, take the 
'thermodynamic' limit — 00 (in a sense of increasing number of zeros — quasi-particles) 
of Gaussian random holomorphic functions (p. The limit — 00 is in fact compatible 
with a semiclassical limit h since the order of a Taylor poljTiomial aproximating 
a Bargmann function (and its zeros) should be at least equal to the number of basis 
states covering the classically accessible area A of phase space (or surface of section), 
N > A/{27rh), while higher coefficients aN+n, n > 0, vanish rapidly with increasing n. 
One can also consider cases of different geometry, for example of kicked spin j systems, 
where quantum (eigen) states can be exactly represented in terms of the so-called SU(2) 
polynomials of finite order 2j . 

It has been shown recently by Hannay for the case of complex coefficients and 
supplemented by Prosen [T^ for the case of real coefficients that the statistics of zeros of 



Gaussian random polynomials are exactly solvable and all k-point correlation functions 



can be given in terms of simple analytical formulae. It has been demonstrated |Tl|] that 
the results obtained by random polynomials indeed reproduce the statistics of zeros of 
Husimi representations of chaotic eigenstates. 



It is the aim of this paper to introduce the local parametric statistics of zeros of Husimi 
representations of the quantum eigenstates. Let us take a family of Hamiltonian systems 
which smoothly depend upon an external parameter A. Then the Bargmann represen- 
tation of a given eigenstate f{z,X) and its zeros Zn{X) are also smooth functions of the 
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parameter A. Therefore we can introduce the velocities f„ as derivatives of zeros Zn with 
respect to an external parameter A 

«„(A) = ^.„(A) = -gj^-j^^ (2) 

Then we define a parametric k-point correlation function of k zeros z = {zi, . . . , Zk) and 
k velocities v = {vi, . . . ,Vk) as 

~p,{z,v) = {6'\z-z')6'\v-v'))^,^^, (3) 

where {)z' v' represents an (ensemble) average over all k-tuples of zeros z' and correspond- 
ing velocities v' of a given state (or ensemble of states). In other words, Pk{z, v)(P'^z (P'^v 
is a probability to find a k-tuple of zeros and corresponding velocities in a small 4fc-cube 
of volume dP^zdP^v around point {z,v). Integrating out the velocities one should obtain 
the usual /c-point correlation function Pk{z) of zeros only |12 



p(z)= Jd''vp{z,v). (4) 

We show below, in section 2, that this parametric correlation functions can be explic- 
itly calculated for Gaussian ensembles of parameter-dependent random polynomials with 
either complex or real coefficients. In section 3, we shall verify our conjecture that the 
obtained result on parametric statistics of roots of statistical ensembles of random poly- 
nomials may be applied to quantum chaotic systems by presenting some numerical results 
obtained in a generic chaotic system, namely 2-dim semi-separable oscillator In 
section 4, we write a closed system of 'equations of motion' for the zero-flow Zn{X) for the 
simplest, linear parametric dependence of the coefficients a„, and stress the integrability 
of the underlying dynamical system. 



2 Parametric statistics of roots of Gaussian random 
polynomials 

2.1 Complex coefficients 

In this section we study parametric statistics of statistical ensembles of Gaussian random 
polynomials. In the first subsection we are dealing with the case of complex coefficients 
while slightly more complicated case of real coefficients will be dealt with in the next 
subsection. Let us write a random polynomial of order in a form 

/(^,A) = f:a„(A)^" (5) 

n=0 

where the coefficients a„(A) (for a fixed realization) may depend smoothly on an external 
parameter A. We will need only first derivatives with respect to A so the parametric 
statistical Gaussian ensemble of random polynomials is completely specified by fixing A 
and saying that a„ and d\an are complex Gaussian random variables with prescribed 
covariances 
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which need not be further specified for the statement of the general result. Fixing the two 
/c-tuples of complex numbers z and i;, we define the 3A; linear combinations of 2(A^ + 1) 
Gaussian random variables a^, d\an-, n = 0, . . . , N 

N 

fj = f{zj,\) = Y.anz], (6) 

n=0 

N-1 

f; = dJiz„X)=Y.in + l)an+izl (7) 

n=0 

d - - ^ 

fj = ttM,A)=/, + V;, fj = dxfiz,,X) = Y.dxa^zJ (8) 

which are again Gaussian random variables. The joint probability distribution of 3A; 
random variables $, = (/, /', /) can be therefore written as 

where M = M(z, v) is 3k x 3k Hermitian covariance matrix which is written in a block 
form as 

/A B D\ 

M = {^(g)C)= C E (10) 

VD^ F / 

= {fjfn, B,, = {f,fr), c,i = {f;fr), (n) 
= (f,f:), E,, = {f'f;), F,i = {hn). (12) 



Note that tilded symbols are used to denote matrices (or vectors or scalars) which explic- 
itly depend on the parametric velocities v. fc-tuple z are the zeros if / = 0, and fc-tuple 
V are the velocities if in addition / = {d/dX)f = 0. So the parametric /c-point correlation 
function (|^) can be written as a linear transformation of a joint distribution P{$) 

p(z,v) = |<iV||^m/',0) (13) 

where d{f, f)/d{z, v) = Y[j=i \fj\^ is the Jacobian of the mapping (z, v) — > (/, /), and 

L = C - B^A'^B - (E - B^A-iD)(F - D^A-^D)-i(E^ - D^A^^B) (15) 

is the central k x k block of the inverse of covariance matrix, M ^. The dependence on 
positions of zeros z and velocities v is digged in the definitions of the matrices (|T0|,|Tl],|T2|) 
In general, the dependence on velocities can be made explicit in the following way. Writing 
a diagonal velocity matrix as V = diag {vj,j = 1 . . .k} and using a definition (^) one can 
observe that the covariance matrices D, E, F have a simple velocity dependence which can 
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be written in terms of their velocity-independent counterparts D, E, F 
D = D + BVt, D,i = {fj;), 

E = E + CVt, E^, = {rj:), (16) 
F = F + Etyt + VE + VCVt, F,i = lfJ*] 



This relations can be used to prove that the determinant of the covariance matrix does 
not depend on velocities 

detM = detM (17) 

where 3fc x 3A; matrix M is obtained from M by replacing the blocks D, E, F by D, E, F. 
Using some elementary algebra one can rewrite the matrix L in the form which makes 
velocity dependence explicit 

L"^ = G"^ + (V^ + G-^K)(H - K^G-^K)-^(V + K^G"^) (18) 
where we have introduced the matrices 



(19) 



G 


= C- 


B^A 


^B 


= Gt, 


H 


= F- 




Id 


= Ht, 


K 


= E- 


B^A- 







The representation of parametric correlations pk in terms of moments of a Gaussian ([14|) 
is very convenient since it may be explicitly evaluated as the sum of all possible pairwise 
contractions of integration variables /j (Wick theorem) and expressed in a compact form 
following an approach of Hannay 

where permanent of a square matrix per S = J2p Hj Sjp. is a symmetric analog of a de- 
terminant det S' = Y,pi.~YT{j Sjp- where p are permutations with signatures (— )^. Inte- 
grating out the velocities, which can be done by putting expression (0) into eq. (|I4D and 
evaluating the inner Gaussian integrals in terms of new variables Uj = vjfj, one obtains 
the fc— point correlation function of Hannay 

Pfc(^) = -j-nnr- (21) 



IT 



det A 



The formula ( pOD is a general result on parametric statistics of Gaussian random poly- 
nomials with complex coefficients. Its important feature is purely algebraic dependence 
on velocities in contradistinction with e.g. parametric energy level statistics (see e.g. 0, 
chapter 6) where velocities have Maxwellian distribution since the Hamiltonian of the 
energy level flow can be clearly written as the sum of the usual kinetic and potential part. 

In the important special case where (the coefficients of) the random polynomial and 
its parametric derivative are statistically uncorrelated 

(andxam) = (22) 
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we obtain parametric correlation functions which are invariant under the change of sign 
of velocities pk{z, —v) = pk{z, v) since E = D = K = and therefore 

L"^ = G"^ + V^P-^V. 

In other words, the average velocity (and all its odd moments) is zero (v) =0. In this 
case, the determinant of 3k x 3k covariance matrix has also a simple factorization in terms 
of k X k matrices 

detM = det AdetGdetF. 

Since the 2— point parametric correlation function depends on 4-complex arguments 
it may be useful to define also the 2-point velocity-moments of the parametric correlation 
functions 

Uvivrvr) = t ^ , / d\,d\2 v'lvivrvrp^izi, Z2, v,, (23) 

which are different from zero only ifk + l = m + n and finite if k + l + m + n< 4. The 
nontrivial velocity-moments, which still depend on the positions of two zeros zi and Z2, 
can be calculated using the Wick theorem from the representation (|1^. Let us quote the 
results for the driftless case (|22|) 

(viv;) = G2iFu/{Gi,G22 + GuG2i), (24) 
Ivil'H') = iFnF22 + F,2F2i)/iGnG22 + GuG2i), (25) 
'vfvf) = 2F^,/{GiiG22 + GuG2i). (26) 



So far the correlations between the coefficients of random polynomials an,dxan have 
been completely arbitrary! Now we specialize to the case where coefficients of random 
polynomials are 5— correlated 

(anflm) = Smnbn, {dx^ndxa*^) = 5mJ)n (27) 

where the variances 6^,6^ are still arbitrary. Introducing two polynomials 

N N _ 

g{s) = Y.bns\ -g{s) = Y.bns'' (28) 

n=0 n=0 

the relevant matrices can be expressed as 

A^i{z) = g{z,z*), (29) 

Bjiiz) = z,g'{z,z!), (30) 

Cji{z) = g'{z,z-n + z,ztg"{z,z:), (31) 

Fji{z) = -giz.zt). (32) 

Putting k = 1, 1— point parametric statistics can be explicitly written in an elegant 
factorized form 

/ N I^iv/(J(Z)) , , 

p^{z,v) = p,{z) ^ (33) 
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where 

TT as as 

is a general density of zeros as can be deduced from | 



and 



is an universal velocity distribution normalized to a unit mean square and cr^(z) 
is a mean square velocity which is inversely proportional to the density of zeros 



(34) 
(35) 

(36) 



So, the theory of random polynomials predicts a universal form of a velocity distribution 
(35) when it is locally rescaled to a unit mean square local velocity. 

In case of eigenstates of RMT in the Bargmann representation one has 6„ = l/n!, 6^ = 
cr^/n! and N — > oo, so 

g{s) = exp(s), g{s) = a^exp(s) (37) 

and therefore the density distribution and the local mean square velocity are constant, 
Pi{z) = l/vr, (|f P) = 0"^, so one has 



In this probably the most important particular case we are also able to give some details 
of 2-point parametric correlation function P2{zi, Z2, f i, V2) which is only a function of the 
4 real quantities instead of 8: distance between roots \z2 — zi\, magnitudes of the velocities 
Ml = |fi|,'U2 = \v2\ and the angle between velocities = arccos(Ref 1^2/1^1^2!) • Writing 
s = 1-22 — -Zip and fixing the velocity scale by putting cr = 1 we may express 2— point 
parametric correlations as 

4(e^-l)V(/?2 + 72^4/37) 



P2 



(39) 



7r2(/5-7)5 
e^' - {s^ + 2)e' + 1 

e'{{e' - l){e' - 1 - s) + aul){{e' - l){e' - 1 - s) + auj) 
e^'ie' - ifie-' - 1 + sf + c?u\u\ - 2(xe\e' - \){e-' - 1 + t)uxU2 cos(0) 

Since this expression is quite complicated it is worthwhile to study its asymptotics for 
large and small distances a/s between the roots. The first two nonzero terms of the small 
s expansion are 

48 



a 
/? 
7 



P2 



-^2 + 





2-8 


Vl 


- V2 




+ \vi + V2 


2)6 



+ 0{s' 



(40) 



TT^{2 + \VI + V2\^)^ 

For large s asymptotics in the leading term, as we expect, the two velocities are uncorre 
lated, while we give also the next term of an expansion in powers of 



P2 = -7 



IT" 



1 + Qe' 



;i+|t;i|2)3(l+|t;2 



2^3 



(41) 



+ 





^-2)(h 


^-2) + 12s{\v,\ 




V2 


2 _ 


- 1 


) - 9S f 1 + V2 


^-9\VI-V2\^ 


(1 + 




2)4 


(1 + 


V2 


2)4 



O e 



-2s 
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Note that the small s expansion pO| ) should be understood strictly pointwise, while it 
is not termwise integrable with respect to velocities vi and f 2 as should be the case for 
the entire 2-point parametric correlation function p2 It is useful to give also the 
velocity-moments (^) which in this case depend only on the distance between zeros i/i 
and probably still contain a lot of information about parametric 2-point statistics 

(v^v;) = -a^e'ie' -l)ie-' -l + s)/uj, (42) 

|^i|>2|') = a\e' + l)ie'-l)Vu, (43) 

'vlvf) = 2a\e^~l)yu;, (44) 

UJ = e'{e' + e^'{e-' ~l + sf. (45) 



2.2 Real coefficients 



In this subsection we discuss the case of parametric statistics of the Gaussian random 
polynomials with real coefficients, i.e. a„ and d\an are real Gaussian random variables 
with prescribed covariances 

) , {andxam) , (dxandxam) ■ 

Fixing the two k-tuples of complex numbers z and v we define 6k real random vari- 
ables Re/j, Im/j, Re/j, Im/j, Re/j, Im/j, or equivalently, a vector of 3 x 2/;; variables ^ = 

(/, /) where f = {f„ f*, ...,/,, /,*), /' = {f[, f[\ . . . , fi /D, / = (/i, A% . . . , fl, h) 
The joint distribution of $, is now again a (real) Gaussian with various blocks and their 
derivations of covariance matrices which are now 2k x 2k matrices and are defined by 
the same formulae ([T0| , p!ll , p!2| , [T6| , [T8| , [T9|) . Using straightforward approach which follows the 
previous subsection and the derivation of nonparametric statistics for real coefficients 
one derives the general formula for the parametric fc— point correlation function of zeros 
of Gaussian random polynomials with real coefficients 



Pk{z,v) 



(27r 



i2fc 



/detL 
detM 



sper 



L2kl 
\L2kl 



L2kl 
L2kl 



Ll2k 
Ll2k 

L2k 2k 
L2k 2k 



Ll2k \ 
Ll2k 



L2k 2k 
L2k 2k 



(46) 



which reduces to the corresponding nonparametric A;— point correlation function 

sper G 



Pk[Z) 



(27r)'=Vdet A 



(47) 



after velocities are integrated out. Semi-permanent of a square 2m x 2m matrix (intro- 
duced in [O]) is a homogeneous polynomial of order m of the matrix elements 



sper S 



E En 



(48) 



jl<-..<jm p r=l 

i, <...<l„l 



where p runs over all m\ permutations of m indices {1 . . . m} and addition of indices jr+fTi 



should be understood modulo 2m. Note that sper 



(g) L ) 7^ sper I L 
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whereas per |^(g)L^= per (g) ^ | 1 ) ) ' ^ matrix on the RHS of 

eq. (^6|) cannot be written hke in the complex case (^). 

Since the case of real coefficients is much more complicated than the case of complex 
coefficients we consider only the final and most important specialization, that is of har- 
monic random polynomials with real coefficients in the 'thermodynamic' limit oo, 
(a„am) = Snm/n\, {andxam) = 0, {dxandxam) = o'^6nm/n\,g{s) = exp{s),g{s) = a'^exp{s). 
Then it is easy to see (like in the nonparametric case |]12[) that the two cases of real 
and complex coefficients are different only when some of the zeros Zj is close to the real 
(symmetry) axis, whereas in the opposite case, where all Zj go away from the real axis, 
the parametric statistics for the case of real coefficients converge to parametric statistics 
for the complex coefficients. Explicit formulae for 1 and 2 point parametric correlation 
functions for this case are too lengthy to write out, we shall only give an asymptotic 
expansion for 1— point function close to the real axis (small y = Imz) 

Mx + ^y,u + ^w) = ^ " ,2^,2 ^^u^/^ + ^ ^^^^ 

It is interesting to calculate the 'root mean square' (RMS) velocity which now depends 
on the distance from the real axis y 

The fact that RMS velocity has a singularity around s = 0, (Ifp) oc s~^, is compatible 
with finding that the small s expansion of the parametric density ( ^91) is not termwise 
integrable with respect to the velocity. 



3 Numerical results 

The predictions of the theory of random polynomials have been verified on the para- 
metric statistics of zeros of Husimi representation of eigenfunctions of a generic chaotic 
autonomous Hamiltonian system with two freedoms which depends on an external pa- 
rameter. We have chosen a recently introduced semi-separable oscillator for which an 
extremely efficient quantization scheme exists based on a construction of an exact quan- 



tum Poincare mapping M, 10|. The Hamiltonian of the system is the following 



H = -\h^dl - \fedl + \{x - asign(y))2 (51) 

where the configuration space is a strip (x,y) G (—00,00) x with hard walls and 

therefore Dirichlet boundary conditions for a wavefunction at y = h^^h^. The system is 
integrable for a = while for a > the potential becomes discontinuous and the system 
becomes increasingly chaotic with increasing a. It has been shown in that the system 
is fully chaotic (with possible islands of stability being negligibly small) for the following 
values of parameters: a = 0.25,6^ = 4, 6| = —11,E = 0.5. For these values and for h = 
0.0003 we have calculated a stretch of 16 consecutive eigenstates with sequential quantum 
number (according to Thomas-Fermi rule) ~ 17 684 000. Then we have calculated a second 
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set of eigenstates according to variation of parameters a = 0.25, 6| = 4 + A, 6| = —11 + A 
with A = 5- 10~^. For each of the eigenstates of the two sets we have calculated the Husimi 
representation on the surface of section {x,px) its zeros (only inside classically allowed 
(chaotic) region of phase space) and calculated the variations of zeros w.r.t. small variation 
of the parameter from A = to A = 5 ■ 10~^. Since the variation of the external parameter 
5A = 5 ■ 10~^ was small enough, the identification of corresponding zeros has always been 
possible and the velocities were numerically well defined. Assuming the conjecture that 
the parametric statistics of zeros of Husimi representation of a chaotic eigenstate should 
be described by the theory of random polynomials, only two parameters remain which 
define the scales, namely the average density of zeros pi (which should be constant inside 
the classically allowed chaotic region of surface of section) and the root mean square 

velocity a = \J (|f p). The density or the number of inside zeros is roughly constant 
for our 16 states while RMS velocity a, which is determined by the sensitivity of a given 
eigenstate with respect to the variation of external parameter A, might exhibit substantial 
fluctuations. Therefore, if we want to merge numerical data of all eigenstates together in 
order to improve statistics, we should first rescale the data for each individual eigenstate 
to a unit RMS velocity. 

In figure 1 we plot such numerical velocity distribution for our stretch of 16 eigenstates 
which contains velocities for 16 x 1800 = 28800 zeros (only the zeros whose distance from 
the boundary of classically allowed region of surface of section was larger than few mean 
spacings were taken into account). It is evident that the agreement with the prediction 
of theory of random polynomials ( p5D is statistically significant. 

In order to compare significantly the numerical estimate of the parametric 2— point 
(or higher) statistics with the prediction of the theory of random polynomials one needs 
a data of a higher statistical quality than those of our numerical experiment. This may 
be easier in some explicit one- dimensional quantum maps and we leave it as a challenge 
for a future work. 



4 Zero-flow as an integrable dynamical system 

In this section we show explicitly that the flow of zeros of a polynomial with respect to 
some external parameter A interpreted as a fictitious 'time' can be written in terms of a 
closed set of 'equations of motion'. We argue that the underlying dynamical system is 
completely integrable and construct explicitly the complete set of integrals of motion as 
functions of zeros and velocities. This construction should be considered as a polynomial 



analog of Dyson- Pechukas- Yukawa p|, |T3[ level dynamics of matrices. 

Let us consider a polynomial of order whose dependence on the external parameter 
A is linear 

N 

fiz,\) = l[iz-z,{\)) = foiz)+\Mz), (52) 

N N 
n=0 n=0 

where without essential loss of generality we have assumed that 

a% = 1, ajv = 0. (53) 
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Writing the velocities Vj = {d/dX)zj = —fi{zj)/ Y\i^j{zj — zi), differentiating with respect 
to A again and using Lagrange interpolation formula for fi{zj) and f[{zj) one arrives to 
the simple closed set of equations for 2N complex (or 4N real) dynamical variables, zeros 
Zj and velocities Vj, j = 1, . . . ,N 



d 



VjVk 



(54) 



From the construction it is obvious that the underlying dynamical system should be 
completely integrable. Writing the polynomial 



N 



N 

E( 

n=0 



/(z,A)= I|(^-z,) = >:(-l)"c„(zx 
in terms of fully symmetric homogeneous functions 

Cn(^l, • • • , Zm) 

one easily sees that their derivatives 

^niZl-, 



zn)z 



N-n 



(55) 



\A\=n 

A(i{i,...,N}3eA 



(56) 



d 



Zn.Vi, 



,Vn 



dX 



zn 



N 



(57) 



are just equal to the coefficients of the perturbation polynomial (see eqs. 

h = (-!)"<_„ (58) 

and therefore independent nontrivial (complex) constants of motion, for n = 1, . . . A^. It 
seems more difficult to give a simple Lagrangian or Hamiltonian formulation of the equa- 
tions (pll). Since the system is completely integrable its Lagrange function is not unique. 
However, it is easy to show that eqs. (0) are completely equivalent to Euler-Lagrange 
equations for any Lagrange function L(zi, . . . , zn, f i, . . . , Z]si) which can be written in 
terms of a symmetric (real) and nondegenerate holomorphic function of constants of mo- 
tion and their complex conjugates 



F(Ji, . . . , /at, , . . . , /^) 
det 



F(/i, . . . , In, I^, . . . , Ilf), 



(59) 



+ 0. 



5 Conclusions 

We have defined a new type of local parametric statistics of quantum eigenstates, which 
describes the statistical distributions of zeros of Husimi functions and the velocities — 
the derivatives of zeros with respect to an external parameter. It has been conjectured 
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that apart from the scahng parameters (density of zeros and mean square velocity) the 
parametric statistics of zeros of a Husimi representation of an eigestate of a classically 
chaotic system should behave universally and should be well described by the ensembles 
of Gaussian random polynomials to the same extent as the statistics of energy spectra, 
eigenvector components, or matrix elements are described by the ensembles of Gaussian 
random matrices. The parametric statistics of roots of Gaussian random polynomials (or 
zeros of Gaussian random holomorphic functions) have been solved analytically and its 
applicability to simple chaotic Hamiltonian systems has been demonstrated numerically. 

In case where the dependence upon external parameter is linear the flow of zeros 
of an analytic (Bargmann) function has been formulated in terms of fully integrable 
dynamical system, where the external parameter and the zeros are playing the roles of 
fictitious time and mutually repelling quasi particles, respectively. The analytical results 
on local parametric statistics presented in the first part of this paper are nothing else but 
'equal-time' correlation functions of such a dynamical system. The fact that the velocity 
distribution (^) is a generalized Lorentzian and not Maxwellian is due to nonseparability 
of kinetic and potential energy. For example, the 2— body case [N = 2) in the center of 
mass frame v = vi = —V2 = 2z, z = Z2 — zi can be described by the Hamiltonian 

H=l\z\'\z\' = M/\z\\ 

where p = \z\^z is a (complex) momentum conjugated to the (complex) coordinate z. 

The generalization to the calculation of the most general 'different-time' correlation 
function — non-local parametric A;— point correlation function 

~p{z',v',\) = (j[6{z'^ - z,{\,mv'^ - {d/dX)z,{\))^ 

is straightforward by putting 

/i = f'j = dzf{zj, Xj), fj = {d/dX)f{zj, \j) 

instead of (P,|^J^) and following the formalism of section 2. 
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Figure captions 

Figure 1: The numerical velocity distribution of the chaotic semiseparable oscillator 
(full curve) (see text for details) as compared to the theoretical velocity distribution (|35|) 
based on the theory of random polynomials (dashed curve). The agreement is highly 
statistically significant. 
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